Non-sequential triple ionization in strong fields 
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We consider the final stage of triple ionization of atoms in a strong linearly polarized laser 
field. We propose that for intensities below the saturation value for triple ionization the process 
is dominated by the simultaneous escape of three electrons from a highly excited intermediate 
complex. We identify within a classical model two pathways to triple ionization, one with a triangular 
configuration of electrons and one with a more linear one. Both are saddles in phase space. A 
stability analysis indicates that the triangular configuration has the larger cross sections and should 
be the dominant one. Trajectory simulations within the dominant symmetry subspace reproduce 
the experimentally observed distribution of ion momenta parallel to the polarization axis. 
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I. INTRODUCTION 

Multiple ionization is a fundamental process in laser- 
atom interactions plQl. It received particular attention 
when it was realized [0 that the cross section was much 
larger than could be expected for independent electrons 
model 11,01. One, therefore, came to distinguish two 
modes of multi-ionization, a sequential ionization pro- 
cesses that is compatible with an independent electron 
model, and a non-sequential one, where correlations in 
electron dynamics are important. The way interactions 
come in is a matter of debate but the rescattering model 
|Q,|| has a fascinating combination of quantum tunneling 
and classical dynamics of electrons in a field and has the 
strongest support from experimental observations, theo- 
retical analyses and numerical simulations [0-EO|. One 
consequence of the rescattering model is that the final 
decay towards multiple ionization is preceeded by the 
formation of a highly excited complex of electrons close 
to the nucleus. This compound states forms during the 
collision of the rescattered electron with the core. In our 
analysis we take this state as the initial condition for the 
final decay towards a multiply ionized atom. 

The way two electrons escape from this compound 
state in non-sequential double ionization has been elu- 
cidated in a recent set of experiments by Weber et al [pj . 
They found that the electrons escape preferentially with 
the same momenta. This observation prompted us to 
analyze the classical pathways that could lead to double 
ionization |gl|,g2[. In an extension of Wanniers analy- 
sis p3|-p7| we found that for two electron escape most 
classical paths lie near a subspace of electrons moving 
with the same distance from the nucleus but reflection 
symmetric with respect to the field axis. Calculations 
of ion momenta parallel and perpendicular to the field 
within this subspace [HJI^ are in favorable agreement 
with experiment. 

At these laser intensities also more than two electrons 
can escape simultaneously. In particular, triple ionization 
has been investigated by Moshammer et al. [0 and the 



recoil ion momentum distributions have been measured. 
It is our aim here to present the extension of the pre- 
vious arguments HJ22] to the case of triple ionization. 
We identify saddles for non-sequential three electron es- 
cape and the relevant symmetry subspaces in which the 
saddles are situated. Numerical simulations within the 
subspaces allow us to obtain the ion momenta that re- 
sult from simultaneous triple ionization. Specifically, in 
section II we discuss the pathways that lead to triple ion- 
ization, identify two saddle configurations and analyze 
their stability properties. In section III we then present 
numerical results for the distributions of ion momenta 
within these two subspaces and compare to the experi- 
mentally observed distributions. Finally, in section IV we 
draw conclusions and given an outlook to higher multiple 
ionization. 



II. PATHWAYS TO NON-SEQUENTIAL TRIPLE 
IONIZATION 

As in the case of double ionization we can assume 
that the electron dynamics during the final decay is fast 
compared to the changes in the field [pl^,p2|. Then the 
straightforward extension of the double ionization anal- 
ysis [ plp2) to three electron escape calls for an investi- 
gation of stationary configurations of three electrons in 
an atom exposed to a strong static field. Keeping a high 
symmetry between the electrons an obvious candidate is 
a configuration where electrons are placed at the vertices 
of an equilateral triangle whose plane is perpendicular to 
the field axis (see Fig. |l|). Deviations from this plane, in 
particular displacements in field direction, will again be 
amplified by electron repulsion, thus leading to electrons 
being pushed towards the nucleus and not to triple ion- 
ization. This configuration has C^^ symmetry and will 
be analyzed in subsection 
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Further investigations show, however, that there is a 
second stationary configuration. This second configura- 
tion has all electrons in a plane with one electron on the 



field axis and the others symmetric with respect to it. 
This configuration i s of symmetry C2v and wiU be ana- 
lyzed in subsection II B , The decay signatures of both 



configurations are different and the experimental obser- 
vations can be used to decide between them. 




Before we proceed with our analysis note a scaling sym- 
metry of the classical Hamiltonian (|l|) that can be used 
to eliminate one parameter. If the variables are rescaled 
according to 

H = F^I^H' 
p = ^i/4p' 



t = F-^/H' 



(4) 



the dynamics becomes independent of the peak field am- 
plitude, i.e. the system is described by the Hamiltonian 
(|l|) with F = 1. We will use this scaling, but drop the 
prime to denote the scaled variables, for the analysis of 
the saddles, but we will keep the full field dependence for 
the dynamical simulations. 



A. C^v symmetry subspace 



FIG. 1. Schematic representation of electron motion in the 
Csv symmetry subspace. The electrons are placed in a plane 
perpendicular to the field polarization axis at the vertices of 
an equilateral triangle with the field axis in the center. The 
circle at the origin of the coordinate system marks the posi- 
tion of the nucleus. 

In both cases we start from the Hamiltonian for three 
electrons in a neutral atom in a laser field polarized along 
the z-axis. Within dipole approximation, with infinitely 
heavy nucleus and in atomic units, the Hamiltonian is 
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and the time dependence of the pulse 

f{t) = s.in^{TTt/Td) cos{ujt -f- 



(2) 



(3) 



F, Td, UJ and stand for peak amplitude, duration, fre- 
quency and phase of the external field, respectively. 

The projection of total angular momentum onto the 
polarization axis is conserved so that the number of de- 
grees of freedom is reduced by one, leaving eight plus the 
time dependence from the external field. We work in a 
subspace where this component vanishes. 



For zero total angular momentum projection on the 
polarization axis, the Cs^ symmetric configuration cor- 
responds in cylindrical coordinates to electron positions 
given by Zi = Z , pi = R and ipi = 27ri/3 with conju- 
gate momenta p^. = Pz/^, Ppi = PrI'^ and p^^ — 0, 
respectively. The Hamiltonian of the system in the Ca^ 
subspace then reads 



H{pR,Pz,R,Z,t) 
with potential energy 
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(5) 



(6) 



The three terms in Eq. (jg) are the cumulative interaction 
of three electrons with the triply-charged nucleus, the re- 
pulsion energy between electrons at distances 2i?sin(7r/3) 
and the interaction with the external field, respectively. 
The variables are scaled according to (^ so that the field 
amplitude F is absent. 

Electron motion close to the nucleus is much faster 
than the change of the phase of the laser field applied in 
the experiments [[7|- pT| , pl 22 1. Therefore, we can use an 
adiabatic approximation and gain insight into the quali- 
tative features of the ionization process by analyzing the 
potential (||) for fixed external field. Note, however, that 
we use the full time dependence for the determination 
of the final ion momenta below. The potential (0) for 
a given time and thus a fixed value of f{t) is shown in 
Fig. g. The saddle is located along the lines Zg = r^ cos 9s 
and R^ — r, sin0, with 9^— 9 ov 9^=1: — 9 where 
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V2 



(7) 



and 



r',=V6/\f{t)\. 
The energy of the saddle is 



(8) 



(9) 



During a field cycle the saddle moves in from infinity 
along the line 9^ = 9, back out to infinity and then in 
and out again along the line 9g = tt — 0. 
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FIG. 2. Equipotential curves of the potential energy (n) 
for three electrons The potential is plotted for a field strength 
Ff{t) — 1. The saddles move along the dashed lines when 
the electric field points in the positive z direction and along 
the one obtained by reflections on z = during the other half 
of the field cycle. 

For the dynamics in full phase space the stability prop- 
erties of the configuration are important. Since the elec- 
trons move fast compared to the frequency of the field 
we can for this analysis employ an adiabatic approxima- 
tion and work with a constant field. Then we can scale 
out the field amplitude completely and put f(t) = I. 
The triangular electron configuration then becomes sta- 
tionary and the dynamics near this saddle is governed 
by the eigenvalues of the linearization. Within the Cs^ 
symmetry subspace the triangular electron configuration 
has one unstable direction corresponding to a crossing 
of the saddle. In these scaled units the Lyapunov expo- 
nent for motion along this reaction coordinate papS] is 
Xr = 1.1054. The other eigenmode of the saddle in the 
subspace is stable. 

Harmonic approximation in the full eight-dimensional 
configuration space results in three additional pairs of 
degenerate eigenvalues. Of these pairs two are stable 
and one is unstable. The unstable eigenspace points to- 
wards asymmetric configurations, has a Lyapunov expo- 
nent Va = 1.4496 and is spanned by the vectors 



P2 = Rs+Wa, Z2 = Zs + 2.0642w;a, 
P3 = Rs- 2.0000?«a, Z3 = Z,- A.l284wa, 



(10) 



and 



Pi = Rs - Ua, zi = Zs - 2.0642Uo, 

P2= Rs+ Ua, Z2 ^ Zs + 2.0642Ua, 

P3 — Rs, Z^ — Zs- 



(11) 



The first mode describes a motion where for positive Wa 
two electrons move away from and one towards the nu- 
cleus or vice versa for negative Wa ■ The other mode leaves 
one electron on the saddle and brings one closer and one 
further out. 

One may compare the period of the laser field applied 
in the experiments [O] with the time scales for crossing 
of the saddle in the C31, subspace and for departure from 
the subspace. For peak field intensity 1.5 • 10^^ W/cm^ 
and for wavelength of 795 nm, the field period (in the 
scaled variables (0)) is 27r/a; = 33.6 while l/\r = 0.9 
and l/i'a — 0.69. This indicates that crossing the saddle 
and departure from the non-sequential triple ionization 
manifold take place rather quickly, justifying the adia- 
batic analysis a posteriori. 



B. C2v subspace 

In addition to the triangular configuration there is one 
with all electrons in a plane. For fixed external field the 
potential (^) possesses a stationary point with one elec- 
tron on the polarization axis and the other two placed on 
either side symmetrically with respect to the field axis 
(Fig. g). Without loss of generality we may assume that 
the electrons are confined to the x-z-plane and then the 
C2V subspace is spanned by 



2/1 = 2/2 = ys = 

Pvi ^ Py2 = Pys = , 
xi = , Pj,^ = , 
X2 = -X3 = X, Px2 
Z2= Z3^ Z, p^2 = 



= -Pxi =Px/2, 

PZ3 = Pzl2 . 



(12) 



The phase space for motion in this subspace is six- 
dimensional, with variables (x, z, Zi,p2,,pz,p2i). The 
Hamiltonian becomes 
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Pi = Rs 



zi= Zs + 2.0642wa, 



(13) 



In contrast to the C^v subspace the saddle cannot be 
given analytically and has to be found numerically. For 
f{t) — 1, the saddle is located at \xs\ = 1.1607, Zs — 
1.1143 and zi.s = 1.4665 and has the potential energy 
Vs = —7.3902. The configuration is elongated along a 



line perpendicular to the field axis, with axial position 
of the central electron further out than that of the outer 
two. Nevertheless, the distance of the central electron to 
the nucleus is smaller than that of the outer electrons. It 
thus experiences a stronger attractive force and the bal- 
ance of force requires that the other two electrons con- 
tribute a component pointing away from the nucleus, i.e. 
that \zs\ < \zi^s\, as indeed observed. 



rotation around the field axis. The neutral mode is con- 
nected with conservation of angular momentum and was 
eliminated in the C^v case by transformation to polar co- 
ordinates; it shows up in the C2v analysis where cartesian 
coordinates are more convenient since one electron would 
come to lie on the singularity of the polar coordinate sys- 
tem. The first unstable mode with Lyapunov exponent 
i'2b = 0.9024 corresponds to deviations along the y-axis, 



■>i 



-1 



-2-L 



yi = U2b 

J/2 = -0.5752M2b 

2/3 = -0.5752M26- 



(16) 



The other unstable mode, with the Lyapunov exponent 
i'2c = 1.3712, is related to symmetry breaking within the 
cc-z-plane. 



XI = 0.0746u2c 

X2 = Xs + 0.6598u2c 

xs = —Xs + 0.6598u2c 



Zl = Zl^s 
Z2 = Zs + U2c 
Z3^ Zs- U2c 



(17) 



The Lyapunov exponents are of similar magnitude as in 
the C31, case and the adiabatic approximation remains 
justified. 



C. Comparison between the subspaces 



FIG. 3. Saddle in the C2v symmetry subspace. The three 
electrons are placed at the stationary point of the potential 
( [L3| ) (for Ff{t) — 1). The arrows indicate the unstable mode 
of the saddle corresponding to the simultaneous escape of all 
electrons. The circle at the origin of the coordinate system 
marks the position of the nucleus. 

Harmonic approximation in the C2V subspace around 
the saddle reveals one stable and two unstable modes. 
One of the unstable directions, with Lyapunov expo- 
nent A2r = 1.0980, corresponds to simultaneous escape 
of three electrons: in the eigenspace 



Zl = Zl,s + W2r 
X = Xs + 0.6183w2r 
Z ^ Zs + 1.1417w2r, 



(14) 



all three components increase for positive W2r- This 
eigenspace thus corresponds to the reaction coordinate 
in the triangular configuration. The other eigenmode 



Zl = Zi^s + U2a 
X = Xg — 0.2313u2a 
Z ^ Zs — 0.3127u2a, 



(15) 



with Lyapunov exponent V2a = 1.7937 is related to dou- 
ble or single ionization, according to the sign of U2a- 

Harmonic approximation in the full space gives in addi- 
tion to the modes within the subspace three stable modes, 
two unstable ones reflecting divergence from the C2V sub- 
space and one neutral mode connected with an overall 



Analysis of the stationary points in phase space shows 
that there are two pathways to non-sequential triple ion- 
ization. The triangular configuration has a potential 
barrier (for scaled field Ff{t) — 1) of —7.6673, slightly 
lower than the one —7.3902 for the saddle in the C2V 
subspace. Therefore, for increasing energy in the com- 
pound complex, triple ionization via the triangular path 
becomes possible first. The triangular state is also less 
unstable than the C2V configuration: it has a largest Lya- 
punov exponent for symmetry breaking perturbations of 
Va = 1.4496, as compared to V2a = 1.7937. For motion 
along the reaction coordinates the Lyapunov exponent 
of the triangular configuration is A^ = 1.1054 and thus 
slightly larger than A2r — 1.0980 for the C2V configura- 
tion. The relative phase space weight of trajectories near 
the symmetry subspaces is determined by a competition 
between motion along the reaction coordinate and ampli- 
fication of deviations from symmetry pq,E^ . Generally, 
the faster the motion along the reaction coordinate and 
the slower the symmetry breaking the larger the phase 
space that is dominated by the saddle. 

There are no quantitative estimates for the phase space 
region influenced by the saddle, except close to a thresh- 
old, where a suitable extension of Wanniers arguments 
can be employed |23|-|25|| . For triple ionization and in the 
presence of several symmetry breaking modes, the gen- 
eralization |26|j27| ] gives a threshold behavior a{/S.E) oc 
(A£')" with exponent 
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(E-' 



A, 



(18) 



where Xj is the Lyapunov exponent of the reaction coor- 
dinate and the Vi are the ones of the symmetry breaking 
modes. Specifically, for our example here we find for the 
triangular configuration 



as = 2.6228 
and for the C2v configuration 

a2 = 3.7043 . 



(19) 



(20) 



This already suggests that the triangular configuration 
will be the more dominant configuration. Besides the 
threshold behavior the two saddles also differ consider- 
ably in the final momentum distributions of ionizing elec- 
trons and this will be analyzed in the next section. 



III. ION MOMENTA DISTRIBUTIONS 

In the present section we show numerical simulations 
of the non-sequential ionization. We restrict our calcula- 
tions to the symmetry subspaces and compare the result- 
ing final ion momentum distributions. Even though the 
symmetry subspaces are of zero phase space weight and 
are furthermore unstable non-sequential triple ionization 
trajectories in full phase space have to pass sufficiently 
close to the saddles and these subspaces so that they will 
show distributions very close to the ones obtained in the 
subspaces pi| , p2| . We can, therefore, also compare our 
data with experimental distributions. In order to allow 
comparison with experiments we now omit the elimina- 
tion of the maximal field amplitude F by the rescaling 

In the experiments of Moshammer et al. ||11[ on triple 
ionization of Ne ultrashort (30 fs) laser pulses at 795 nm 
wavelength and with peak intensities of 1.5 PW/cm^ were 
used. This corresponds to a frequency of w = 0.057 a.u. 
(atomic units), a pulse duration, measured as full width 
at half maximum, of 11 periods 2tt/u! and a maximal 
field strength of F = 0.207 a.u.. Among the measured 
data we focus on the distributions of ion momenta par- 
allel to the polarization axis. In the limit of negligible 
small momentum transfer by the absorbed photons, the 
ion momentum pi^n reflects the sum of the momenta of 
the emitted electrons, pion = ~^iPi [3'II3- 

The three-electron Hamiltonian (ll) corresponds to 
lithium. In order to relate the calculations to other atoms 
with more electrons interactions with core electrons have 
to be neglected and the point of reference in energy has 
to be shifted to the threshold for triple ionization. Specif- 
ically, for the modeling of the experiments [Q on triple 
ionization in Ne we thus assume that in a rescattering 
process the energy transfer is less than the threshold for 
triple ionization (about 4.6 a.u.). The precise value of 
the energy transfered in the rescattering event and thus 
of the energy E of the initial compound state cannot be 
determined within our model and constitutes a free pa- 
rameter. Similarly, the time io during the pulse when the 



rescattering complex forms is a second free parameter. 
But both parameters can be determined rather reliably 
by comparison with experiments 123]. We now turn to 
the simulations within the two symmetry subspaces. 
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FIG. 4. Final distributions of parallel ion momenta cal- 
culated in the Czv symmetry subspace for the initial en- 
ergy E — —0.5 a.u., peak field amplitude F = 0.207 a.u., 
and pulse duration T^j = 20 x 2tt/u!, where lo — 0.057 a.u.. 
The initial times under the pulse envelope are to = 0.25rd 
(a), to = 0.75Td (b), to = OAT^ (c), and to = O.GTd (d). 
Dashed lines indicate the estimates ±3F/p(to)/aj. Note that 
the distributions are essentially the same independently if one 
chooses io before or after the peak field value provided /p(to) 



is the same, 
trajectories. 



The distributions are based on about 8 • 10 



A. Simulations in the C3„ subspace 

In a first series of simulations we focus on the triangu- 
lar configuration. We fix the initial energy E and analyze 
the dependence of the results on the initial time ip . To 
this end we chose a microcanonical distribution of ini- 
tial conditions in the C^y symmetry subspace for energy 
E = —0.5 a.u. and random phases (j) for different initial 
time to. The final distributions of ion momenta parallel 
to the polarization axis are shown in Fig. H. As in our 
previous analysis of the non-sequential double ionization 
process [ plp^ the distributions has a double hump struc- 
ture and the width of the distributions increases when the 
initial time Iq approaches the peak amplitude of the field. 
The width of the distributions can be estimated from the 
maximal energy a free electron can acquire in the field, 
i.e. twice the ponderomotive energy. For three electrons 
ejected in the same direction, the corresponding maximal 
ion momentum depends on the amplitude of the field at 
the point in time when they are ionized, i.e. 



3F/p(i) 



(21) 



where fp{t) = s\r? {-nt/Td) is the pulse envelop. If this 
time is taken to equal the starting time to of the simula- 
tions we find values of Pmax which correspond very well 
to the widths of the distributions in Fig. ^. Fig. |4| indi- 
cates also that the distributions are basically the same 
independent of whether io is chosen before or after the 
maximum of the pulse, provided fp{to) is the same. This 
implies that the dominant ionization takes place during 
the first field cycle after formation of the complex. Fit- 
ting the width of the calculated distribution to the ex- 
perimental results allows one to estimate the moment in 
the pulse when the majority of the triply ionized ions are 
created. 
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FIG. 5. Ion momentum distributions in the Csv subspace 
for fixed initial time io = 0.5Td and different initial energies 
E = -0.05 a.u. (a), E = -1.5 a.u. (b), E = -2 a.u. (c), and 
E = —3 a.u. (d). The parameters are the same as in Fig. ^ 
The minimal energy of the saddle, obtained for the maximal 
field amplitude, is —3.49 a.u.. 

The dependence of the ion momentum distributions 
on the initial energy E for fixed initial time to is pre- 
sented in Fig. ^. The widths of the distributions does 
not change significantly with E but their shape does. 
For initial energy close to the minimal saddle energy, 
which for F = 0.207 a.u. is —3.49 a.u., the distributions 
show a single maximum only while for higher energies 
two humps form. The electrons that cross the saddle 
when the energy is near the saddle energy are slow and 
the combined interactions of the external and Coulomb 
fields shapes the distribution |22]. The situation is dif- 
ferent for the high energy case: then shortly after the 
electrons cross the saddle the interaction with the elec- 
tric field is stronger than the attraction to the nucleus 
and the distribution is mostly shaped by the laser field 
p2| . As the initial energy of the complex is determined 
in a rescattering event and is higher for stronger exter- 
nal fields, the distribution of ion momenta should show 
a transition from a distribution with a single maximum 
near the threshold for triple ionization to one with a dou- 



ble maximum higher up. For even higher fields the con- 
straints from the triangular configuration are relaxed and 
less symmetric modes of triple ionization become possi- 
ble. 



B. Numerical simulation in the C2v subspace 

We now turn to the ion momentum distributions in 
the C2v subspace. Much of the analysis proceeds as 
in the previous subsection, except for the difference in 
Hamiltonians (here it is (|l^)) and the choice of initial 
conditions. The microcanonical ensemble of the initial 
conditions is not straightforward to realize because (even 
for fixed time) for nonzero external field the system is 
open. In the previous section, the difficulties could be 
overcome by constraining initial conditions to lie on the 
energy shell and in the hypersurface Z = (for details 
see |2^). Here such a restriction is not sufficient, phase 
space remains open. The troublesome configurations are 
those where electron distances to the nucleus are strongly 
asymmetric, one being close and two far away or vice 
versa. Since we assume that the initial state is formed 
during a rescattering and should be confined to region 
close to the core, such configurations cannot be formed. 
We, therefore, choose initial conditions microcanonically 
and require in addition that distances of the electrons to 
the nucleus are not larger than the minimal distance of 
the saddle at maximal field to the nucleus. 



0.08 




FIG. 6. Final distributions of parallel ion momenta calcu- 
lated for non-sequential triple ionization in the C2v symmetry 
subspace for the peak field amplitude F — 0.207 a.u., pulse 
duration T^ = 20 x 27r/aj (where ui = 0.057 a.u.) and ini- 
tial time in the pulse duration to = 0.5Td. The initial energy 
is E = -0.05 a.u. (a), E = -0.5 a.u. (b), S = -2 a.u. 
(c) — the minimal energy of the saddle is —3.36 a.u.. The 
distributions are obtained from a few thousands trajectories. 



As in the previous section, the width of the momentum 
distributions determines the initial time io when triple 
ionization takes place. The dependence on initial energy 
E is also similar, with a collapse of all structures into 
a single peak near zero momentum as the initial energy 
decreases towards the energy of the saddle (Fig. H). 
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FIG. 7. Examples of trajectories corresponding to simul- 
taneous three electron escape by crossing the saddle in the 
C2v symmetry subspace. The electrons, when crossing the 
saddle, move initially towards the positive z direction (the 
saddle in both cases is located at z > 0). However, depending 
on the phase of the field, their final momenta can point in 
either the same (a) or opposite (b) direction. 

In Fig. we show final parallel ion momentum distri- 
butions corresponding to non-sequential triple ionization 
in the C2v subspace for three different initial energy E. 
The structure of the distributions is significantly different 
from that in the previous section. In Fig. p|(a) and 0(b), 
besides two local maxima on the edges of distributions, 
a strong contribution zero momentum value that is split 
into two maxima appears. The origin of these structures 
can be related to qualitatively different electron paths. 
Initially, after crossing the saddle all three electrons es- 
cape towards the same side of the nucleus. However, the 
final momenta along the polarization axis may point to 
either the same or opposite directions, depending on the 
initial conditions and the phase of the field, as shown in 
Fig. M. If the final momenta point in the same direction, 
they can be much larger than if they point in opposite 
directions. Therefore, the edge maxima are the effect of 
the electrons escape towards the same direction while the 
central structure corresponds to an escape in opposite di- 
rections. 



C. Comparison with experiment 



Given the qualitatively very different distributions of 
final ion momenta in the two subspaces we can use a com- 
parison with experimental data |11[] to draw conclusions 
about the dominant ionization mode. In Fig. we show 
the experimental parallel ion momentum distribution to- 
gether with the one calculated from the C^v symmetry 
subspace. In the calculation we have chosen the initial 
time to = 0.33Td that gives the width of the resulting dis- 
tribution in agreement with the experiment. The initial 
energy has been taken as -E = — 1 a.u. but it can be taken 
in a wide range of values above —1.5 a.u. without signif- 
icant changes in momentum distributions (Fig. ph. The 
main features of the experimental distribution are well 
reproduced in our calculations. However, the minimum 
in the distribution is less pronounced than that observed 
in the experiment. Overall, we take this good agreement 
as strong indication that triple ionization occurs in the 
neighborhood of the symmetric process discussed here. 
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FIG. 8. Comparison between experimental (a) and numer- 
ical (b) ion momenta. Ion momentum distributions are from 
the experiment of llll] on triple ionization of Ne atoms in the 
focus of 795 nm, 30 fs (i.e. about 11 x 2tt/u>) laser pulses 
at peak intensity of 1.5 PW/cm^ (i.e. peak field strength 
F — 0.207 a.u.). The numerical distributions are for the Cs^ 
symmetry subspace with initial energy E — —1 a.u. and 
io = 0.33Td where Td/2 = 11 x 2n/uj, see Eq. (||). 

The distribution for perpendicular ion momenta, which 
have been measured in the experiment pd| , can not be 
calculated within the symmetry subspace since this com- 
ponent vanishes exactly by symmetry. In full phase space 
ionizing electrons pass close to the saddle but not exactly 



symmetrically and consequently the transverse ion mo- 
menta do not vanish. 



IV. CONCLUSION 



number of stationary configurations. The identification 
of a dominant pathway will then require an appropriate 
stability analysis and an application of the generalized 
Wannier threshold law |26lj27f|. 



We have considered the pathways along which a highly 
excited three electron complex can decay towards triple 
ionization. In non-sequential multiple ionization of atoms 
in strong laser fields such a complex forms when the 
rescattered electron collides with the core. The subse- 
quent non-sequential decay then has to follow the path- 
ways discussed here. We found two possible pathways, 
one proceeding via a symmetric triangular configuration, 
and one via a planar elongated configuration. A stabil- 
ity analysis shows that the ionization close to the C^y 
symmetry subspace occurs at lower field intensities and 
has a larger cross section (close to threshold). Trajec- 
tory simulations within this subspace give the final ion 
momentum distributions in very good agreement with 
the recent experimental data [Q . The C2v configuration 
has a higher threshold and smaller cross section and gives 
a momentum distribution different from the experimen- 
tally observed one, so that we conclude that the triangu- 
lar configuration is indeed the dominant one. 

We have not discussed pathways to sequential triple 
ionization, and one can expect that there are several, 
including sequential ionization of one electron after the 
other, a double ionization followed by a single ionization 
or the other way around, or various combinations of par- 
tial ionizations and rescatterings. As in the case of double 
ionization [El 22 1 one can also have sequential ionization 



from failed attempts to triple ionize along one of the sym- 
metric saddles. The experimental results of [|l^ suggest 
that, at least in the intensity range considered here, cor- 
related non-sequential electron escape dominates. This 
is further supported by our previous analysis of double 
ionization |23|, where we found that the double hump 
structure observed in the experimental ion momentum 
distributions ||7|-pl| does not appear in the sequential 
process. That is, the electron momenta are not corre- 
lated and the resulting distributions show a strong single 
maximum instead of the double hump structure. 

Finally, let us mention that part of our analysis can 
be easily extended to non-sequential multiple ionization 
where the number of escaping electrons, N, is greater 
than three [£1[. The triangular configuration extends to 
a placement of electrons on a ring with C'nv in plane per- 
pendicular to the field axis. Calculation of the positions 
and the energy of the saddle shows that the energy is non- 
monotonic and begins to increase with N for more than 9 
electrons. For more than 13 electrons the saddle does not 
exist any more. Therefore, for many electrons the config- 
urations with highest symmetry will disappear and others 
will dominate the ionization process. For triple ionization 
we have identified one, the planar configuration, and for 
more than three electrons we expect a rapidly increasing 
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